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Abstract 


The  problem  of  a  steady  boundary  layer  or  sheath  between 
a  plasma  and  a  magnetic  field  is  considered.   A  self-consistent 
transition  layer  is  found  which  joins  a  uniform  magnetic  field 
at  plus  infinity  with  a  collisionless  field-free  plasma  region 
with  arbitrary  velocity  distribution  at  minus  infinity;  i.e., 
a  magnetic  field  profile  is  found  such  that  the  exact  particle 
orbits  in  this  field  produce  a  current  which  gives  rise  to 
this  field.   An  interesting  feature  of  the  solution  is  that, 
with  any  nonsingular  velocity  distribution  at  minus  infinity, 
the  magnetic  field  penetrates  only  a  finite  distance  into  the 
plasma,  although  the  plasma  extends  to  infinity,  exponentially 
attenuated,  into  the  magnetic  field  region.   The  scale  of 
length  is  the  Larmor  radius.   Electric  fields  arising  from 
charge  separation  in  the  case  of  particles  of  different  mass 
are  ignored  in  this  paper. 
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THE  BOUNDARY  LAYER  BETlAfEEN  A  PLASMA  AND  A  MAGNETIC  FIELD  -  I 


1.   Introduction 

An  idealization  which  has  proved  to  be  extremely  useful 
io  that  of  an  Interface  or  current  aheet  separatins  a  perfectly 
conductins  fluid  from  u  magnetic  field.  Equilibrium  Is 
characterized  by  equality  of  the  raasnetic  pressure,  B  /2p.  , 
and  the  fluid  pressure,  p.   Analysis  of  this  type  of  equilibrium 
configuration  can  be  accomplished  by  classical,  powerful  mathe- 
matical tools  ;  analysis  of  stability  was  at  once  the  earliest 
and  most  completely  developed'";  and  in  the  form  of  the  Cusped 

Geometry,  such  a  model  offered  the  earliest  example  of  an 

3 
absolutely  stable  confined  conducting  fluid'^. 

Nevertheless,  there  arc  many  pvirposes  for  which  this  model 
is  inadequate  and  some  measure  of  the  thickness  of  the  transi- 
tion layer  or  even  its  fine  structure  is  required.   This  is 
necessary,  for  example,  to  compute  particle  losses  from  a  Cusped 

This  vrcrk  was  done  largely  in  1955  In  connection  with  the 
problem  of  estimating  particle  losses  in  cusped  geometries. 

1.  H.  Grud,  WASH-184,  p.  l44  (Jan.  1955);  also  H.  Grad  and  H. 
Rub i n ,  Pr'.'C.  Second  International  Conf.  on  the  Peaceful  Uses  of 
Atomic  Energy,  Vol.  '^1 ,   p.  190  (195<^J  ^nd  J.  Berkowitz  et  al. 
Ibid.  p.  171. 

2.  M.  D.  Kru.skal  and  M.  SchwartzschlJd .  Proc.  Roy.  See.  k22Z> , 
548  (195^);  ... 

H.  Grad,  WASH-289,  p.  115  (June  1955)  and  TID-7505,  p.  2^8 
(Feb.  I95S);  J'  Berkowitz  et  al. ,  Proc.  Second  International 
Conf.  on  the  Peaceful  Uses  of  Atomic  Energy,  Vol"  31,  p.  177 

(1958). 

3.  H.  Grad,  WASH-289,  p.  115  (June  1955). 
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Geometry,  to  refine  stability  criteria,  and  to  estimate  electron 
cyclotron  radiation. 

For  simplicity,  we  consider  an  essentially  one-dimensional 
model  in  which  all  quantities  depend  on  the  single  coordinate  x 
(traversing  the  sheath),  and  the  magnetic  field  is  uni-dlrectlonal, 
say  in  the  direction  z,  B  (x).   The  existence  of  an  unambiguous 
solution  of  this  type  requires  several  Idealizations  and 
specializations.   First,  to  allow  such  a  time-independent  solution 
requires  the  neglect  of  dlsslpative  processes  (macroscopically  -- 
resistivity;  microscopically  --  particle  collisions).   Next  we 
observe  that  the  problem,  as  posed,  does  not  have  a  unique 
solution.   In  the  macroscopic  version,  any  pressure  profile 
p(x)   and  magnetic  field  profile  B  (x)   satisfying  the  simple 
condition  p  +  B  /2|j,  =  constant  defines  an  allowable  equilibrium 
solution.   Similarly,  we  shall  see  that  there  is  a  correspondingly 
very  large  variety  of  microscopic  solutions  depending  roughly  on 
how  many  particles  are  placed  on  each  magnetic  line. 

This  ambiguity  can  be  removed  by  asking  for  the  thinnest 
transition  zone  compatible  with  certain  given  properties  at 
infinity.  More  precisely,  it  is  shown  that  there  exists  a 
unique  solution  corresponding  to  any  given  velocity  distribution 
at  minus  infinity,  provided  that  there  are  no  "trapped"  particles 
which  do  not  have  access  to  infinity.   Such  a  solution  can  be 
Interpreted  as  the  microscopic  solution  which  best  represents 
the  originally  described  macroscopically  sharp  interface.   But 
more  important  than  this  purely  formal  property  is  the 
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Glgnlflcance  of  the  thinnest  layer  In  yielding  the  smallest 
possible  particle  loss  rate  in  a  Cusped  Geometry  and  also  the 
loviest  possible  radiation  loss  rate  in  any  contained  plasma. 
It  is  therefore  a  theoretical  Ideal  which  one  should  try  to 
approximate  experimentally. 

It  is  very  easy  to  estimate  the  order  of  magnitude  of 
the  thickness  of  this  layer.   The  idealized  current  sheet 
requires  an  infinite  current  density.   But,  for  a  given  mean 
thermal  speed  C,   the  maximum  current  density  that  can  exist 
in  the  plasma  is  approximately 

J  =  neC. 

Thus,  in  a  layer  of  thickness  5,  the  maximum  change  in  B  is 
given  by 


B  =  l-L  J5  =  M,  neC5, 


Ignoring  constants  of  order  unity,  we  have 

,2 


p  =  nkT  =  nmC 


from  which 


5  ^   B     mC  B^ 


li^neC   eB  M-^p 


Consequently,  a  change  in  B  from  the  value  zero  to  a  value 
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giving  equality  of  magnetic  and  gas  pressure  can  occur  only 
over  a  distance  at  least  as  large  as  the  representative  Larmor 
radius,  mC/eB»  This  argument,  which  makes  no  direct  reference 
to  particle  orbits,  is  consistent  with  the  concept  of  a  Larmor 
radius  as  the  minimum  distance  required  to  turn  back  a  particle 
incident  on  a  magnetic  field.   Furthermore,  we  shall  see  that 
the  statement  of  self-consistency,  viz.,  that  the  computed 
particle  orbits  give  rise  to  the  precise  current  required  to 
produce  the  field  gradient.  Is  exactly  equivalent  to  a  more 
refined  version  of  a  pressure  balance o 
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2.   Particle  Orbits 

We  consider  a  magnetic  field  with  a  single  component 
in  the  z-dlrectlon  and  represent  the  component  unambiguously 


y 


'1 


f*    J 


^A 


B 


V 


FIG.  I 


by  the  scalar  symbol  B.   B  is  assumed  to  vary  in  the 
x-direction  only.   Thus  both  the  vector  potential  and  the 
current  can  be  taken  in  the  y-dlrection,  and  we  represent 
them  by  the  scalar  symbols  A  and  J  respectively;  we  have 


(2.1) 


f          R  -  ^A 

J          1    dB 

1    d^A 
^^0    dx^ 

The  z-component  of  velocity  of  any  particle  is  a 
constant;  we  are  therefore  free  to  consider  the  motion  to 
be  in  the  plane  (x,y).   The  two  dimensional  Hamlltonian 
governing  the  motion  of  a  particle  of  mass  m  and  charge 
e   (signed)  is 


8 


(2.2) 


H  =  liT  tPx  +  (Py  -  ^A)'l 


where 


(2.3) 


Px  = 


P.r  = 


L 


dx 
mu  =  m  ^ 


mv  +  eA  =  m  ^jI^  +  eA 
at 


Since  the  y-coordinate  Is  Ignorable,  the  problem  can 
be  treated  as  one -dimensional  with  the  same  Hamlltonlan  (2.2), 
but  with  p   now  Interpreted  as  a  given  constant  of  the 
motion.   As  a  one-dimensional  system,  the  energy  integral, 
H  =  constant,  immediately  gives  a  general  solution  explicit 
in  terns  of  a  quadrature  involving  the  presumed  given  function 
A(x)  as  well  as  the  two  constant  values  of  H  and  p  . 
However,  almost  all  features  of  interest  can  be  seen  from  the 
energy  equation  itself  before  quadrature;  viz.. 


(2.4) 


P  mu  +  f (x)  =  e 


where 


(2.5) 


IM    =  ^  [l^  -  eA(x)]2. 


We  use  the  symbols  e  and  \i     for  the  constant  values  taken 
by  H  and  p   respectively. 
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The  motion  is  that  of  a  particle  under  the  influence  of 
the  potential  |.  For  example,  if  the  range  of  values  taken 
Toy  A  is  sufficiently  great,  the  inequality  f  <  e  implies 
that  the  entire  motion  lies  within  the  bounded  interval  on  the 
X-axis  determined  by 


(2.6)  [|i  -  eA(x)]^  <  2m€ 


On  the  other  hand,  if  A  remains  bounded  as  x  recedes  to 
plus  or  minus  infinity,  then  for  appropriate  values  of  the 
constants  €  and  \i,   infinity  may  be  accessible  to  the  particle 
in  its  motion.   Such  statements  can  be  made  precise  if  certain 
qualitative  features  of  A(x)  are  known. 

To  proceed  further,  we  restrict  our  attention  to  special 
magnetic  field  profiles  which  are  qualitatively  similar  to  the 
expected  solution  of  our  problem.   Specifically,  B(x)   is 
positive  and  approaches  the  limit  zero  at  minus  infinity  and 
the  finite  limit  B-,^  at  plus  infinity;  in  particular,   J  =  0 
at  both  plus  and  minus  infinity.   We  also  expect  B  to  be 
monotone,  but  this  restriction  is  not  required  to  carry  out  the 
analysis.   If  we  further  assume  (to  be  verified  afterwards) 

X 

that  /  B(x)dx  converges  at  the  lower  limit,  then  A(x)   con- 

—  00 

verges  to  a  finite  value  at  minus  infinity,  which  value  we  can 
take  to  be  zero.   A(x)  increases  monotonely  and  is  unbounded  at 
plus  infinity.   Concisely, 
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(2.7)  • 


A(-co)    =  0, 

A{+oo)      =     +«>, 

A'(x)   >   0 

B(-co)    =   0, 

B{-h^)    =  B^, 

B(x)   >   0 

j(-oo)     =    0, 

j{-Ho)    =    0 

*^ 


It  should  be  emphasized  that  these  are  properties  which 
we  hope  the  solution  will  have.   Whether  they  are  in  the 
nature  of  boundary  conditions  that  can  be  imposed  at  will, 
or  whether  they  may  be  automatically  satisfied  or  are  even 
incompatible  remains  to  be  seen. 

Adopting  these  restrictions  on  A(x),   we  can  now  make 
much  stronger  statements  about  the  possible  orbits.   For  the 
potential  J(x),  we  have  the  three  possibilities  shovm  in 
Figure  2  depending  on  whether  \x     has  the  same  or  opposite 


♦-x 


(b)/ie  <  0 


FIG.  2 


>  X 


(c)>ie  >  0 
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sign  as  e.   In  the  two  cases  |ie  <  0,   ^  is  monotone.   For 
any  value  of  e   consistent  with  the  existence  of  a  real  orbit, 
the  particle  "rolls  up"  the  Incline,  turns  around,  and  returns 
to  minus  Infinity.   The  velocity,   u,   decreases  monotonely 
(In  time)  through  the  value  zero  and  ultimately  approaches  the 
negative  of  Its  Initial  value.  The  turning  point,  u  =  0,  Is 
given  by  the  unique  solution  of  the  equality  (2.6).   In  the 
case  lie  >  0,  a  sufficiently  small  value  of  e  yields  a 
periodic  orbit  in  which  the  two  extremes  of  x  are  given  by  the 
equality  (2.6),   On  the  other  hand, a  large  value  of  €  yields 
an  orbit  similar  to  the  previous  ones,  except  that  u  is  no 
longer  monotone.   We  call  the  periodic  orbits  trapped  and  the 
ones  accessible  to  infinity  free. 

Quantitatively,  the  condition  that  an  orbit  be  free  is 
simply 

(2.8)  e  >  iJ.V2m. 

In  the  case  \ie   <  0,   there  are  no  orbits  which  violate  this 
inequality;  in  the  case  |j.e  >  0,  the  trapped  orbits  are 
characterized  by 

(2.9)  0  <  e  <  |j.^/2m. 

We  shall  find  that  the  particle  orbit  properties  just  given 
are  the  only  ones  that  enter  into  a  full  determination  of  the 
self-consistent  solution.   For  completeness  we  show  in  Fig.  3 
a  representative  set  of  orbits  in  the  (x,y)-plane  and,  in  Fig, 
4,  typical  phase  plane  orbits;  in  the  latter,  \i     is  constant 
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In  each  figure  while  e  Is  varied. 


y^ 


-»  X 


>x 


(c)/ie  >0,  «  >/iV2m 


(d)/ie  >0  ,€<  fji^/Zm 


FIG.  3 


u  A 


(a)  ae  <  0 


(  b ) ^e  >  0 


FIG.  4 
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3.   General  Self-Consistent  Formulation 

We  are  looking  for  a  configuration  or  assemblage  of 
charged  particles  which  will  be  described  by  a  time -independent 
distribution  function  f(x,u,v).   The  dependence  on  the  third 
component  of  velocity,  w,  has  been  omitted  .   We  normalize 
f  to  be  a  n"umber  density  and  note  the  macroscopic  connections. 


(3»1)         n(x)  =  /  f  du  dv 
(3.2)         J(x)  =  e  /  V  f  du  dv 


(3.3)         ^xx^^^  =  m  /  u^f  du  dv. 


The  quantity  P    is  the  Indicated  component  of  the  stress 
tensor.   If  f  is  isotropic,  it  is  exactly  the  pressure,  p. 
It  should  be  noted  that,  in  virtue  of  symmetries  that  will 
appear,  the  macroscopic  mean  flow  velocity  in  the  x-directlon 
is  zero,   /  uf  du  dv  =  0. 

The  distribution  function  f  satisfies  Liouville's 
equation  (sometimes  called  the  collisionless  Boltzmann  equation 
or  the  Vlasov  equation)  which  states  merely  that  f  is  a 
constant  along  every  curve  in  phase  space  which  is  the 
possible  orbit  of  a  particle.   In  this  problem,  it  suffices  to 


4,   To  be  more  realistic,  one  can  consider  f  to  be  an  even 
function  of  w;  this  condition  for  all  x  is  compatible 
with  the  equations  of  motion.  Thus  there  will  be  no  resul- 
tant gas  velocity  or  current  in  the  z-direction,  and  all 
the  results  hold  as  given. 
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take   (x,u,v)  as  the  phase  space  instead  of  (x,y,u,v)  which 
Is,  a  priori,  correct,  since  it  is  easily  seen  that  there  is 
a  unique  orbit  which  passes  through  each  finite  point  (x,u,v).-^ 
This  entire  family  of  orbits  is,  in  principle,  known  if  A(x)  is 
supposed  to  be  given.   A  solution  of  Liouvllle's  equation  will 
then  be  specified  by  assigning  values  to  f  on  each  orbit. 

Now,  we  knov/  that  e  and  ja  are  constant  on  an  orbit. 
This  suggests  that  we  use  them  as  parameters  to  specify  the 
various  orbits.  The  statement  that  f  is  constant  on  an  orbit 
could  then  be  interpreted  as  a  statement  that  f  is  a  function 
of  e  and  la.  This  is  approximately  true,  but  requires  elabora- 
tion.  It  is  not  necessarily  true  that  a  given  pair  of  values 
{e  ,[i)      singles  out  a  unique  orbit.   If  A(x)   is  not  monotone, 
there  xiill  be  distinct  sets  of  orbits  with  the  same  values  of 
e  and  \x;   clearly,  one  may  wish  to  assign  different  values  to 
f  on  the  different  sets  of  orbits.   It  is  easy  to  verify,  how- 
ever, that  if  A(x)   is  monotone,  this  assignment  is  unique. 
Geometrically,  this  is  clear  from  Pig.  2  and  the  related 
discussion. 

To  effect  the  Integration  over  f,  (3.l)-(3.3),  it  is 
convenient  to  Introduce  e  and  ix  as  integration  variables 
instead  of  u  and  v.   From  the  defining  relations 

^     1/2    2 , 
e  =  p-  m(u  +  V  ) 


(3.4; 


2 
M,  =  mv  +  eA 


5.   In  this  context,  the  point  u  =  v  =  0,  x  =  constant  is  to  be 

interpreted  as  an  orbit.   These  singular  orbits  can  be  ignored 
in  the  problem  at  hand  since  they  do  not  contribute  to  the 
current . 
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in  v;hich  A  enters  as  a  parameter,  we  see  that  the  quantities 
n,  J,  and  P   will  be  given  as  functions  of  A  alone.   The 
self-consistency  relation,  (2.1),  takes  the  form  of  a  second 
order  differential  equation,  viz.,  d  A/dx   is  a  specified 
function  of  A.   There  is  an  alternative  procedure,  viz.,  to 
use  the  exact  macroscopic  equation  of  momentiom  conservation, 

(5.5)         P^^  +  2^-  ^   =  constant. 

o 

This  is  a  first  order  equation  expressing  B  =  dA/dx  as  a 
function  of  A  and  is,  or  course,  merely  the  first  integral 
of  the  previous  second  order  equation.   However,  it  will  be 
iliunlnating  to  carry  through  both  procedures  explicitly. 

Before  carrying  out  this  plan,  we  must  specify  the  distri- 
bution function,   f.  This  is  to  a  large  extent  arbitrary. 
For  those  values  of  e  and  \i     which  correspond  to  free  orbits 
extending  to  minus  infinity,  this  specification  amounts  to  a 
boundary  condition  on  f  at  minus  infinity.  Thus  a  complete 
specification  of  f  amounts  to  the  imposition  of  a  boundary 
condition  plus  the  specification  of  the  density  of  trapped  par- 
ticles for  each  appropriate  €  and  p..   The  value  of  |i  locates 
the  potential  minimum,  Fig.  2c,  and  thus  the  approximate  position 
of  the  trapped  particle.   The  separation  into  free  and  trapped 
reglonii  is  shown  in  the   (€,M.)-plane  in  Pig.  5  and  in  the 
(u,v)-plane  for  a  specific  value  of  A  in  Fig.  6.   Fig.  5 
is  dravm  for  eA  >  0|  the  labels  "trapped"  and  "no  orbits" 
are  to  be  Interchanged  if  eA  <  0.   Similarly,  in  Fig.  6,  the 
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parabola   opens  downward  if     eA  <  0. 


trapped 


>   € 


>U 


FIG.  6 


We  now  return  to  the  transformation  (3.^)  from  the 
variables  {u,v)  to  the  variables  (e,M.).   It  is  evident  that 
each  of  the  half  planes  u  >  0  and  u  <  0  is  mapped  into  the 
same  domain  of  (c^M-).   Specifically,  each  half  plane  is 
mapped  one-to-one  into  the  parabola 


(3.6) 


2m6  >  (la-eA)' 


This  multiple  mapping  merely  introduces  a  factor  two  in  the 
Jacobian;  (the  integrand  in  all  the  expressions  {J>.1)-{J>.J>) 
is  even  in  u).   We  compute  (for  u  >  0) 


(3.7) 


mu 

mv 

|(ijj4 

0 

m 

2  2 

=  m  u  =  m[2m€  -  (ti-eA)  ] 


1/2 


Writing  F(€,M.)  for  f(x,u,v)  when  expressed  as  a  function 


-  17  - 


of  the  Indicated  variables  we  obtain 

(3.8)  n  =  I  /  [2me  -(n-eA)'']     F(€,|i)ded|i 

Pp  p  -1/2 

(3.9)  J  =  ^  /  (n-eA)[2me-(|i-eA)'^]    P(€,M.)dedix 

m 

P  P  1/2 

(3.10)  Pxx  =  "T  ^  [2m€-(n-eA)'']   P(e,n)d€dM. 

m 

The  range  of  Integration  is  in  all  cases  given  by  (2.6).  The 
dependence  of  n,  J,   and  P^,^  on  A  enters  explicitly  In 
the  integrand  and  In  the  limits  of  integration.   The  integra- 
tion domain  is  indicated  by  the  dashed  line  in  Fig.  5  (for 
eA  >  0);  it  overlaps  both  free  and  trapped  contributions. 

It  is  important  to  note  that  a  boundary  condition  on  f 

at  minus  infinity  restricts  F(e,|i)   only  within  the  parabola 

2 
2m€  >  [I   .      For  example,  taking  f  to  be  Isotropic  at  minus 

infinity  Implies  that  F  is  a  function  of  e  alone,  not 

everyv/here,  but  only  within  this  parabola.   In  particular, 

the  condition  that  there  exist  no  trapped  particles   (F  =  0 

for  2me  <  |J.  )  precludes  the  possibility  that  F  be  isotropic 

anywhere  but  at  minus  infinity. 

In  (5.10),  P^„  is  expressed  directly  as  a  function  of 

A.   We  ncv/  wish  to  compute  dP^^/dA,   remarking  that  the 

domain  of  integration  as  well  as  the  integrand  depend  on  A. 

However,  since  the  integrand  happens  to  vanish  on  the  boundary 

of  the  domain  of  integration,  the  only  contribution  comes  from 
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differentiation  of  the  integrand,  viz. 


(5.11)    -gf^  =  %   ^  [2me-((x-eA)^]    F(e,n)d€dn 


m 


=  J 


J 


Thus  we  verify  the  equivalence  of  the  two  self-consistency 
formulations 


(5.12)         _i_4!A^j(A) 


^^o  d? 


2 
1   /dA' 


(5-1^)         Pxx(A)  +2^(h|)   =^0 

o 

The  differential  equation  for  the  determination  of  A(x) 
is  exactly  that  of  the  motion  of  a  particle  in  a  given  poten- 
tial j  we  interpret  x  as  time,  A  as  position,  B  =  dA/dx 
as  velocity,  and  M-^P^yC^)  ^^  ^^®  potential  which  is  a  given 
function  of  position.  Thus  we  can  determine  the  properties 
of  the  solution  A(x)  from  monotonicity  and  similar  qualita- 
tive properties  of  the  function  Pyv(A).   What  we  shall  do  is 
discuss  functions  P-^y(A)  which  yield  satisfactory  solutions 
A(x),  and  then  discuss  the  nature  of  the  functions  F(€,m,) 
which  give  rise  to  such  P  „(A). 


6.   Since  we  shall  wish  to  consider  discontinuous  functions 
F(€,M,)  and  the  integral  is  improper,  this  differentia- 
tion formula  should  be  verified  directly;  this  offers  no 
difficulty. 
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First  we  remark  that  the  constant   P   is  the  value 

o 

taken  by  P    at  minus  infinity  (i.e.,  at  A  =  0);  if  f 

happens  to  be  isotropic  at  mlnum  Infinity,  then  F   is 

o 

the  gas  pressure.   The  condition  J  =  0  at  minus  infinity 
is  a  condition  on  F,  viz.. 


/  ix[2m€  -  n  ] 


-1/2 


F(e,|i)ded|j.  =  0. 


This  is  satisfied  if,  for  example,  F  is  an  even  function 
of  \x      (at  minus  infinity,   f  is  an  even  function  of  v). 
More  specially,  this  is  satisfied  if  f  is  isotropic  (P  is 
independent  of  |i)  at  minus  infinity.   For  simplicity,  we 
shall  henceforth  consider  F  to  be  even  in  ij,;  this  is  no 
restriction  on  trapped  orbits  since  they  occur  for  only  one 
sign  of  [I. 

We  have  already  determined  that  P„„(0)  =  P   and 
p'  (0)  =  J(0)  =  0;  (the  accent  denotes  a  derivative  with 
respect  to  A).   We  also  wish  to  have  P„„(oo)  =  0;  (there  is 
no  gas  at  plus  Infinity).   The  simplest  possibility  would  be 
to  have  P„^(A)  monotone  (Fig. 
7).   In  this  case  the  solution  P^x^ 
to  the  differential  equation 
(3.15)  evidently  has  a  mono- 
tone analogue  "velocity",       _ 
B  =  dA/dx,  which  increases 
from  the  value  zero  at  A  =  0 
to  the  limiting  value  B-j_  given 


FIG.  7 


■>A 
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by 

(3.14)        B^/2|x^  =  P^ 

as  A  approaches  infinity.  Now  consider  x  to  decrease  from 
plus  infinity;  i.e.,  the  "ball  rolls  up  the  hill"  in  Pig.  7. 
As  A  decreases  to  zero,  both  B  =  dA/dx  and  J  =  -(d  A/dx  )/|i, 
approach  zero.   In  the  "normal"  case  that  the  curvature  of  the 
potential  function  is  finite   [P„„  (0)  is  finite],  then  Just 
as  for  a  pendulxam,  it  takes  an  infinite  amount  of  "time"  to 
reach  the  summit;  A  =  0  is  attained  at  x  =  -00.   On  the  other 
hand,  singular  behavior  in  the  form  of  an  infinite  curvature 
can  yield  the  unusual  result  A  =  dA/dx  =  d  A/dx  =0  at  a 
finite  position,  x  .   In  this  case,  the  complete  solution  to 
the  problem  consists  in  piecing  together 


(3.15) 


A  =  0,  X  <  x^ 


A  >  0,   X  >  Xq, 


Since  there  is  no  magnetic  field  to  the  left  of  x  =  x  ,  the 
particle  orbits  are  rectilinear  and  all  gas  properties  are 
constant.   In  other  words,  it  is  irrelevant  whether  we  apply 
the  boundary  condition  at  x  =  -00  or  at  x  =  x.  or  at  any 
intermediate  point.  To  confirm  that  this  analytically  artifi- 
cial function  solves  our  problem,  it  is  only  necessary  to  note 
that  the  orbits  computed  in  this  proposed  field  are  smooth 
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(the  acceleration  is  at  least  continuously  differentiable)  and 
do  Indeed  yield  the  appropriate,  self-consistent,  current. 

To  complete  the  discussion  it  is  necessary  to  find  what 
conditions  on  F  guarantee  that  P^„(A)  is  monotone  and  under 
what  conditions  P''  (O)  is  finite  or  infinite.  This  is  done  by 
the  following  three  theorems  which  are  proved  in  Appendices  I 
and  II. 

Theorem.  1.   If  F(€,|i.)  is  even  in  p.,  and  for  each  e,  F 
is  a  monotone  decreasing  function  of   ||i|,  then  P^x^^^  "^^ 
monotone  decreasing. 

In  the  special  case  of  an  isotropic  f  at  minus  infinity, 
this  requires  that  the  number  of  trapped  particles  of  a  given 
energy  decreases  as  we  move  to  the  right.  More  specially,   P^^ 
is  monotone  for  an  isotropic  distribution  at  infinity  and  no 
trapped  particles. 

Theorem  2.   If  F(€,|J,)   is  a  sufficiently  smooth  function  (e.g., 
continuously  differentiable),  then  J  is  of  order  A  for 
small  A  and  P"  (O)  is  finite. 

Theorem  J>,      If  F(€,m,)   is  discontinuous  across  the  curve 
separating  trapped  from  free  particles,  then  J  is  of  order 
/aT"  for  small  A. 

In  particular,  an  isotropic  distribution  at  infinity  with 
no  trapped  particles  is  necessarily  discontinuous. 
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4.   calculation  of  the  "Thinnest"  Profile 

We  now  consider  the  special  case  in  which  there  are  no 
trapped  particles. 


(^.1) 


F(€,n)  =  0,   e  <  nV2m. 


The  simplest  way  to  incorporate  this  condition  is  to  cut  off 
the  integration  at  e  =  p,  /2ra.   Combining  this  with  the 
inequality  (3-6),  introducing  the  new  variables  a  and  6 
(as  in  Appendix  I), 


r 


(4.2) 


a  =  2m€ 
6  =  |i  -  eA 


and  taking  e  >  C  (there  Is  no  loss  of  generality),  we  obtain 
the  range 


(^.3) 


n 


-a  <  6  <  a-eA 


a  >  -^  eA. 


We  have 
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a-eA 


(4.4) 


a-eA 


(a^-  C^)-^/2  ^(^^  6+eA)d6 


(a^-  d^)"-'-/^  0{a,    6+eA)dC 


(a^-  C^)^'^^  jZ^(a,  6+eA)d6  (  , 


where  we  have  set  j6{a,[i)     for  F(a  /2m,  |j,). 

For  the  special  case  of  an  isotropic  distribution  at  minus 
infinity,  j2J  is  a  function  of  a  alone,  and  one  integration 
can  be  performed  explicitly. 


(^.5) 


■xx=  "T  I  j(a-eA)(2aeA-eV)    + 


eA 


2r7r 


-1 


+  a  [J  +  sin"-^(l-eA/a)]>  jZ^(a)ada. 


In  this  case,  the  monotonlcity  with  respect  to  A  (cf.  Appendix 
I)  can  be  Inferred  directly  from  the  monotonlcity  of  the 
integrand. 

From  Appendix  II,  we  take  the  behavior  for  small  A, 


(4.6) 


J  =  -XA-*-/^  +  0(A) 
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where 


00 

r 


(4.7) 


X  =  -% 


ni  J, 


{2ea)^^'^   j2J(a,a)da   ; 


-jZJ(a^a)  Is  the  Jump  in  0{a,\i.)  as  described  in  the  Appendix. 
For  an  isotropic  distribution,  iZ^(a,a)  is  merely  0ia,).  From 
this  we  conclude. 


:4.8) 


p   -  p  = 

XX     o 


»  ^xP^^  +   O(A^). 


The  solution  of  the  differential  equation  (5,13)  in  the  neigh- 
borhood of  A  =  0,   viz 


^(«)'   2xaV2, 


2[i      ^dx 
^o 


is  then 


(4.9) 


r 

A 

<     B 

J 


We  see  that  not  only  the  current  but  even  its  derivative  is 
continuous  at  the  join. 

Possibly  the  most  interesting  special  case  is  that  of  an 
isotropic  Maxwellian  distribution  at  minus  infinity. 
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mn. 


(4.10) 


< 


P(^)  =^iW-  ^^p(-^At^) 


mn. 


^("^  =  2B^  expC-aVSnikT^). 


The  expression  for  P    is  not  appreciably  simpler  except  that 
the  scaling  can  be  shown  explicitly. 


4   r~ 
(4.11)    P^/p^  =  ^s)  =  1^  J  QO)  exp(-  I  pV)pdp 


1/2 


where 


1/2 


2rTr  .  _.„-l, 


(4.12)    Q(p)  =  0-l)(2p-l)    +  p^[5  +  sln--^(l  -  l/p)] 


and 


(4.15)    s  =  eA/(mkTQ) 


1/2 


The  dimenslonless  vector  potential  s  can  be  Interpreted  as 
distance  measured  in  local  Larmor  radii. 


(4.14) 


S  =  eB/(mkT^)'  ' 


The  differential  equation  for  s  takes  the  form 


1  /ds 


(4.15)        ^  (^)  =1  -Pis) 


where 
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r 


i  =  x/A 


(4.16) 


a2  = 


m 


•^o  o 


T  • 


Except  for  scaling,  there  is  essentially  one  magnetic  profile 
to  be  computed.  The  scale  length,  "K,      can  be  interpreted 
variously  as  the  ratio  of  the  speed  of  light  to  the  appropriate 
(ion  or  electron)  plasma  frequency,  or  as  a  Larmor  radius  com- 
posed of  the  magnetic  field  strength  at  plus  infinity  and  the 
mean  thermal  speed  at  minus  infinity. 

The  constant  X.  ,  can  be  evaluated  explicitly  (cf.  (4.7) 
and  (4.10)),  and  we  find  that  near  x^j 

(4.17)    B/B;^  =  -^  r   (5A)(^-^^)  ^   [0,125(^-^0)]   ' 

917" 

the  initial  scale  length  is  about  eight  times  A. 

Certain  singular  distributions  exhibit  interesting  mathe- 
matical features.   For  example,  a  two-dimensional  constant 
speed  distribution. 


(4.18) 


mn. 


F(e)  =  ^  5(e-€Q) 


allows  P    to  be  computed  explicitly. 


(4.19) 


Pxx=  0   ' 


2  eA  >  a^ 


Pxx  =  ?  ^o^o((l-^A/°^o)(2eA/ao   -   e^A^/a^) 


1/2 


+  5  +  sin^"(l  -  eA/a   )  >     ,         i 


2  ^  "•""      \-   -  -"/^oM     '  ■?  —    ^  -o 


eA  <  a^ 
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It  Is  worth  remarking  that  any  truncated  distribution  function, 
i.e.,  one  which  is  zero  for  sufficiently  large  energy,  exhibits 
this  feature,  that  P    vanishes  for  sufficiently  large  A. 
Such  a  distribution  will  therefore  yield  a  finite  transition; 
both  plasma  and  field  are  together  in  only  a  finite  layer. 

A  more  interesting  singular  distribution  is  one  which  is 
singular  in  \i, 


(4.20) 


0{a,\i)   =  V'(a)5(n). 


This  contains  particles  all  of  which  are  initially  aimed 
perpendicular  to  the  magnetic  field. 

An  elementary  computation  yields  the  exact  results 


(4.21) 


^  n=  2 


r^ 


m  Ji 


2   2  2  ~  > 
(a  -  e  A  )    •|{'(a)ada 


'eA 


J  =  - 


2 
2e^A 


m' 


"  -1/2 

2   _2,2.  ^' "^ 


J 
eA 


(a  -  e'^A  )    ^(a)ada  =  - 


efA 
m 


n. 


and  the  expansion 


(4.22) 


xx 


=  P. 


^a2  + 

2m  ^  +  •• 


Thus,  for  this  very  singular  velocity  distribution,  P  is 
regular  near  A  =  0  and  the  magnetic  field  extends  to  minus 
infinity,  decaying  exponentially.   If  T^(a)  is  the  one- 
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dimensional  Maxwelllan  distribution. 


(4.25)        f   =  n^m^^ 


1/2 


exp(-a  /2mkTQ) 


then  we  have,  explicitly. 


(4.24) 


n  =  n  exp(-e  A  /arnkT^) 

J  =  -(e^n^A/m)  exp(-e^A^2mkTQ) 


P^x  =  Pq  exp(-e^AV2inkTQ). 


This  result  for  a  distribution  which  is  singular  in  the 
y-component  of  velocity,  v,  is  extremely  special.   If, 
instead  of  a  concentrated  distribution  at  v  =  0,  we  take  a 
distribution  concentrated  at  v  =  v  ?^  0  (to  maintain  J  =  0 
at  infinity  we  take  the  sura  of  two  such  concentrated  distribu- 
tions at  V  =  ±  V  ),  then  we  return  to  the  general  result, 

l/p 
J --^  A  '  ,  and  the  magnetic  field  no  longer  penetrates  to  minus 

infinity. 
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Appendix  Is  Monotonicity  of  P„^(A) 

We  assume  that  P(e,p.)  Is  an  even  function  of  |j,  and  Is 
a  non-lncreaslng  function  of   I|j,|.  We  wish  to  show  that  J  Is 
negative  for  A  >  0.  This,  of  course.  Implies  that  P    Is 
monotone.  The  expression  (5.9)  for  J  can  be  written 


(I.l) 


J  = 


2e 

m-^  J 


0        f  /5A+a 
ada  < 


eA-a 


2       2  -^/2 
(M,-eA)[a  -(pL-eA)  ]    ;3(a,M.)d|x( 


where  we  have  Introduced  the  variable  a  Instead  of  €, 

(1.2)  a^  =  2m€, 

0(a,M.)  for  F(a  /2m, ij,),  and  the  explicit  range  of  Integration 
(3.6).   Setting 

(1.3)  d  =  M.  -  eA, 

the  Inner  Integral  becomes 


9-  _i/2 

6(a^-  d^)    /2J(a,eA+6)d6 

-a 

r"    .    .  -1/2 


6(a^-  6^)     [0(a,eA+C)-0(a,eA-5)]dd 


The  bracket  In  this  last  expression  has  a  sign  opposite  to  that 
of  eA.   For  eA  >  0,   |eA+d|  >  |eA-d|,   thus  iZf(a,eA+d)  < 
j^(a,eA-d),  and  conversely  for  eA  <  0.   Noting  that  there  Is  a 
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factor  e  in  the  expression  for  J,  we  conclude  that  J  Is 
non-positive.  This  establishes  that  P    is  monotone  decreasing 
in  the  weak  sense.  For  a  given  value  of  A  >  0,  the  inequality 
will  hold,   J  <  0,  unless  0{a,eA+6)   =  0{a,eA-6)      for  all  a  >  0 
and  all  6,  0  <  e  <  a.   In  other  words,  J(A)  =  0  only  if  0 
is  a  function  of  a  alone  in  the  range  eA-a  <  p.  <  eA+a,  i.e., 
f  is  isotropic  for  all  particles,  free  or  trapped,  whose  orbits 
are  accessible  to  the  given  value  of  A.   If  the  distribution 
of  free  particles  is  non-isotropic,  then  J  <  0  for  every 
A  >  0  and  there  is  strict  monotonicity  of  P^^.   Or,  if  the 
distribution  of  free  particles  is  isotropic  but  there  is  a  strict 
inequality  between  the  number  of  trapped  and  free  particles  for 
any  value  of  the  energy,  again  there  is  strict  monotonicity  of 

^xx* 
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Appendix  II:  Behavior  of  P^^(A)   for  Small  A 

Introducing  the  same  variables  as  In  Appendix  I,  we  have 


2e 


(II. 1)    J(A)  =^   ada  I  J  6(a^-  6^)  iZJ(a,eA+6)d6[ . 

"^  o     /  -a 


If  0     is  a  sufficiently  regular  function  of  its  second  argument 
then  so  is  J(A)|  successive  derivatives  can  be  obtained  by- 
differentiation  under  the  integral  sign.   In  particular,  if 
jS     is   differentiable  with  respect  to  p,  then  the  derivative 
j'(0)  =  p"   (0)   exists. 

Now  let  us  consider  the  case  of  a  function  0{^,[i)     which 
is  discontinuous,  specifically  across  the  parabola  e  =  M-  /2m 
which  separates  free  from  trapped  particles.  For  simplicity 
we  write 


0  =  0Q  +  0^ 


where  0       is  regular  over  some  large  domain  and  0-.      vanishes 
on  the  free  orbits,  €  >  n  /2m.  We  have  two  contributions  to  J. 


J  =  J^  +  J^, 


where,  as  before,  J  is  at  least  of  order  A,  For  J-,  we 
evaluate  only  the  trapped  orbits  (for  simplicity  taking  the 
case  e  >  O), 
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J-.  = 


2e 


ado' 


r^ 


2   .2 


-1/2 


S(a  -  6  )    i2^-,^(a,eA+6)d6 


a-eA 


For  continuous  J^-,,  the  inner  integral  can  be  approximated 
near  A  =  0  by 


-1/2 


6(a^-  6^)    dd  =  ^-^(a,a)(2aeA-e^A^) 

1/2 
--*  i3^(a,a)(2aeA) 


1/2 


a-eA 


1/2 
This  gives  the  dominant  contribution  to  Z ,     of  order  A  '  , 


o 


The  quantity  jZl^(a,a)   is  merely  the  jump  in  )Z^(a,|i)   on  the 

1/2 
parabola  p.  =  (2me)  '    ,   or -straight  line  p.  =  a. 
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